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ABSTRACT 

Electron acceleration in high Mach number perpendicular shocks is investi- 
gated through two-dimensional electrostatic particle-in-cell (PIC) simulation. We 
simulate the shock foot region by modeling particles that consist of three compo- 
nents such as incident protons and electrons and reflected protons in the initial 
state which satisfies the Buneman instability condition. In contrast to previous 
one-dimensional simulations in which strong surfing acceleration is realized, we 
find that surfing acceleration does not occur in two-dimensional simulation. This 
is because excited electrostatic potentials have a two-dimensional structure that 
makes electron trapping impossible. Thus, the surfing acceleration does not work 
either in itself or as an injection mechanism for the diffusive shock acceleration. 
We briefly discuss implications of the present results on the electron heating and 
acceleration by shocks in supernova remnants. 

Subject headings: supernova remnants - shock waves - plasmas - cosmic rays - 
acceleration of particle 



1. Introduction 

The discovery of synchrotron X-rays from young supernova remnants (SNRs) provides 
the evidence that el ectrons are accelerated to highly relativistic energy in SNR shocks 



( iKoyama et al.lll995l ). Since SNR shocks are coUisionless, various processes of particle ac- 



celeration may tak e place. A l though the most popular a cceleration mechanism is diffusive 



shock acceleration (lBelllll978l : iBlandford fc Ostrikerlll978l ). electrons are not expected to be 



accelerated easily because they behave adiabatically in the field made by ions in the first 
order approximation. Thus, in order for electrons to enter the acceleration process they 
should be pre-accelerated to relativistic energy by some injection mechanisms. Surfing ac- 
celeration at perpendicular shocks has been studied and regarded as one such prospective 
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mechanism, as well as a n efficient direct acceleration mechanism ( iHoshino fc Shimadal 12002 
McClements et al.lboOll ). 



The basic mechanism of electron surfing acceleration for high Alfven Mach number 
perpendicular shocks is as follows. The shock refiects some of the inciden t ions to th e 
upstream, where the foot region forms on a spatial scale of ion gyroradius (ILeroyl Il983l ). 
The plasma in the f oot region consists of incident ions and electrons and refiected ions. 
Papadopoulod (119881 ) proposed that in this situation, incident electro ns and refiected ions 
excite electrostatic waves by the Buneman instability ( iBunemanl Il958l ) because the relative 
velocity between them is fairly large compa red with the electron thernial ve locity and because 
the waves heat electrons non-adiabatically. ICargill fc Papadopoulod (Il988l ) performed a one- 
dimensional hybrid simulation and demonstrated that strong electron heating actually does 



occur. 



To examine the electron momentum distribution, IShimada fc Hoshind (120001 ) performed 
one-dimensional full-PIC simulations in which both electromagnetic and electrostatic waves 
were taken into account but with a small proton/electron mass ratio. Their results showed 
that significant electron surfing acceleration occurs, which can be explained in the following 
way: in the shock rest frame, electrostatic waves excited by the Buneman instability prop- 
agate away from the shock front with a large amplitude in order to trap some fraction of 
electrons that are accelerated in the direction of Vsh x S by the motional electric field (Vsh 
is the shock velocity and B is the magnetic field). Since the potential structure is uniform 
in the direction of this acceleration in one-dimensional simulations, electrons are kept to be 
trapped and accelerated indefinitely. The degree of acceleration is determined by the coher- 
ence length of the potential structure to the acceleration direction and the central question 
is whether it is long enough for real shocks in SNRs. 

In this Letter, we address this question by using a two-dimensional electrostatic PIC 
simulation. In §2 we describe the setting and results, followed by a discussion in §3. 



2. Simulation 

To perform two-dimensional simulations, we confine our attention to the foot region 
through a proper modeling instead of solving the whole shock structure. Our simulation box 
is taken to be at rest in the upstream frame of reference, i.e., that of incident protons and 
electrons. We do not solve electromagnetic waves; we concentrate on electrostatic waves and 
adopt a real proton/electron mass ratio of 1836. 
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2.1. setting 

We define the a:-direction as the shock normal pointing to the shock front, and thus 
the refiected protons move in the — x-direction. The magnetic field is taken to be spatially 
homogeneous pointing in the 2;-direction and we solve the particle motion and electric field 
in the x — y plane. As the initial condition, we prepare upstream electrons, upstream 
protons and reflected protons. Each population is uniformly distributed in the x — y plane 
and their momentum distribution is given by a Maxwellian at the same temperatures T — 
Te = Tp = Tr = 7eV. In addition, reflected protons have an extra drift velocity in the 
x-direction of = —0.04c (V^ = 2V^h)- Number densities of each population are taken 
as Ue = 1.25np = brir = lcm~^, where subscripts e, p and r represent upstream electrons, 
upstream protons and reflected protons, respectively. These parameters are typical of young 
SNRs. 

We employ the periodic boundary condition in both the x- and y-directions. The electric 
field is solved by a Poisson equation. We have examined various values from to 72 fiG for 
magnetic field strength (uJce/^pe = ~ 0.08), where LUce = eB/nieC and Upe = {Aimee^ /rrieY/'^ 
are electron cyclotron frequency and electron plasma frequency, respectively. Among them, 
in this Letter, we present only one specific case for which the maximum energy of accelerated 
electrons is largest in the one-dimensional (ID) simulation; i.e., the magnetic fleld strength 
is 27//G (u;ce/'^pe=0.03). 

The size of the simulation box to the x- and ^/-directions is taken to be = Ly = 16ABun, 
with a total of 256 x 256 cells, where Asun = '^'^^pe/Vd is the wavelength of the most unstable 
mode of the Buneman instability. Thus, the length of each cell Ax— Ay is 3 times the Debye 
length. We also perform ID simulations for the purpose of comparison with two-dimensional 
(2D) simulations. The number of macroparticles is taken so that initially each cell includes 
10,240 and 80 electrons in the ID and 2D cases, respectively. Time step At is taken as 
O.Ola;"^ and the simulation is followed until lOOOu;"^. 



2.2. behavior of the electric field 

First, we discuss the time development of the electric fleld. The energy density of the 
electric fleld is shown in Figure 1. The three curves show the spatially averaged energy 
density of the electric fleld. The solid curve shows the ID simulation, while the dashed and 
dotted curves show the x- and y-components in the 2D simulation, respectively. These three 
curves show a comparable growth rate although the growth in the 2D case appears to start 
earlier, probably due to a smaller number of macroparticles per cell used and larger initial 
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fluctuations. It reaches maximum at about a hundred times WpJ-, with an e- folding time 
of about lOcjpg^, as the hnear theory predicts. It should be noted that the amplitude of 
the ^/-component is comparable to that of the x-component, which is discussed below. The 
electric field strength for 2D simulation is a few times smaller than that of the ID case, and 
its time history is far smoother. This is because the phase of electron plasma oscillation is 
less coherent in the 2D simulations than that in the ID case. Two-dimensionality of electric 
field structure is clearly seen in a landscape representation of the electrostatic potential 
at if: = 150c<j~g^, as shown in Figure 2. The typical wavelength in the x-direction is Abuh, 
while that in the y-direction shows a larger dispersion and the typical coherence scale of the 
potential fluctuations in the y-direction Lcoh turns out to be on the order of Abuid too. 

The fact that obliquely propagating modes have a c omparable growth rat e with parallel 
propagating modes was pointed out as early as 1960 by iBludman et all (119601 ). who derived 
the dispersion relation against obliq uely propagating ele ctrostatic two-stream instability for 
cold beams without a magnetic field. iLampe et al.l (119741 ) discussed thermal effects and made 
numerical simulations to study nonlinear evolution for initially cold unmagnetized plasma 
and a somewhat smaller ion-electron mass ratio. Because the magnetic field is sufficiently 
weak in our simulations, we may safely ignore effects of the magnetic field on the growth 
timescale of the most unstable mode. The dispersion relation for the present case, ignoring 
thermal effects, is written as 







UJ 



pr 



where tUpr = (47rnre^/mp)^/^ and we safely ignored the effect of upstream protons and the 
relativistic effect. The wavenumber in the x-direction of the most unstable mode is given by 



27r 



A 



(2) 



Bun 



Note that equation (1) is the same as the dis persion relation derived by Buneman for a ID 
case and does not contain ky (lBunemanlll958l ) . This means that if k^ satisfies the instability 
condition, th e system is un s table against all values of ky for cold beams, as previous studies 
have shown. iLampe et al.l (1l974l ) showed that as electrons are rapidly heated, waves with 



ky comparable with Ag^^ survive. Our 2D simulation results are fully consistent with their 
results; growth of the y-component of the electric field is a natural feature of the two-stream 
electrostatic instability. We note that this natural feature is not taken into account in ID 
simulations that discuss the surfing acceleration. 
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2.3. momentum distribution of electrons 

Momentum distribution of electrons in the 2D simulations is also different from that in 
the f D case. Figure 3 shows the spatially averaged momentum distribution of electrons at 
t = 690a;pg^ We find that in the ID case some electrons are accelerated to the —Vy direction 
(and subsequently rotate counterclockwise), as seen in several jetlike features in the velocity 
space; these are identified with the surfing acceleration. Therefore, we confirm that surfing 
acceleration occurs in the ID simulation as has been discussed in the literature. In contrast, 
in the 2D case, we do not see any jetlike structures as seen in the ID case. The distribution 
is concentrated in a central round region and thus no surfing acceleration occurs. 

Figure 4 shows the energy spectrum of electrons at the end of the simulations. A high- 
energy component representing the surfing acceleration is seen in the ID case, while no high 

energy components are seen in the 2D case. The energy distribution is smooth but is not 
described as Maxwellian with a single temperature. In Figure 4, the distribution for the 2D 
simulation without a magnetic field is also shown for comparison. The result with a magnetic 
field is somewhat hotter than that with no magnetic field. A question may be raised as to 
whether the heating increases with magnetic field strength and whether we get sufficient 
acceleration for a sufficiently strong magnetic field. To see this, we estimate the effect of the 
magnetic field on electron energy change over a coherence scale of potential Lcoh, as 

AE = e-y^i^coh = 27rameVi (—] , (3) 

C \UJpeJ 

where we set 

Lcoh = aAeun, (4) 
with a being of the order of 0(1). This is an order of magnitude too small to explain the 
difference in the simulations because ujce/^pe is 0.03. When we inspect the time history of 
electron distribution, we observe that the difference appears after the electric-field strength 
begins to decay. So we think that the difference is caused by the thermalization process 
after the saturation of electric field fluctuations. Actually, our simulation results for vari- 
ous magnetic fleld strengths do not show a signiflcant difference in the energy distribution. 
Therefore, we conclude that surflng acceleration does not work in the 2D simulations or in 
real cases. 



2.4. interpretations 

In essence, in reahstic 2D cases, the coherence length of waves excited by the Buneman 
instability is the order of Abuh in both the x- and y-directions, with a negligible contribution 
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of very long wavelength modes required for surfi ng acceleration. Our sira ulations do not 
show stochastic surfing acceleration mentioned by iMcClements et al.l (120011 ) . either. This is 
because the stabilization of finite amplitude waves by BGK mode does not work in the 2D 
case; there are no large-amplitude waves, since electrons leave the potential on the timescale 
of the electron cyclotron period. Even when they are trapped in the x-direction, they leave 
it to the ?/-direction, which hampers nonlinear growth of the waves. 

Thus, the net result is not the acceleration but the heating, as was discussed in the 
original investigations of Papadopoulos (1988). The electron temperature in the final state 
in our 2D simulation , inferred frorn the o btained velocity dispersion, is about Tf ~ meVj, the 
same as the result of iLampe et al.l (119741 ) without a magnetic field. Proton temperature does 
not change much during the simulation because protons react on a much longer timescale. 
This temperature corresponds to the threshold temperature for the Buneman instability. As 
was shown in §2.3, electric-field amplitude decreases after it attains the maximum. This does 
not mean the electric field energy is simply transformed into the electron thermal energy. 
The former is much smaller than the latter. The electric field plays a catalytic role in 
transforming the kinetic energy of reflected protons. As the electron temperature increases, 
the growth of the electric field is suppressed, and when it attains the threshold temperature, 
the wave growth ceases. The detailed e xaminations o f the heating process are interesting 
but beyond the scope of this Letter (see iLampe et al.l (119741 )). 



3. Discussion 

The present study has made some assumptions that we beheve are not so crucial as 
to invalidate the physics conclusions. First, we have limited ourselves to only electrostatic 
modes and do not solve electromagnetic waves. This is sufficient in the linear regime but may 
influence nonlinear evolution of the instabilities. Although the two-stream situation is known 
to be also unstable against electromagnetic modes, their growth rate is orders of magnitude 
smaller than that of electrostatic modes. In the late stage when the electrostatic modes 
saturate and begin to decay, electromagnetic modes may play an important role in electron 
heating and acceleration. The behavior will surely depend on the electron temperature 
attained in the initial stage dominated by the electrostatic modes. 

In this Letter, we have presented the result of only one specific parameter set. The 
results for other parameter sets will be discussed elsewhere. Here, we briefly discuss effects 
of the initial electron temperature. We have discussed the case in which the drift veloc- 
ity is about 7 times the thermal velocity of the electrons (T = 7eV). We have performed 
simulations for other values of the initial electron temperature, keeping the drift velocity 
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unchanged. When the temperature is 4 times lower, i.e., T — 1.75eV, the structure of poten- 
tial in the ^/-direction becomes more fluctuated because modes with shorter wavelengths in 
the ^/-direction can grow, consistent with our argument in §2.2. Conversely, when the initial 
temperature is higher, growth of short- wavelength modes is suppressed, and long-wavelength 
modes preferentially survive; thus, one may expect that more heating or even surfing accel- 
eration may be realized. In reality, this is not the case. We have done simulations for which 
the temperature is 4 and 16 times higher, i.e., T = 28 and 112 eV. The resultant coherence 
length of potential in the y-direction becomes longer, as expected, but it is only a few times 
longer and does not much change the situation. It should also be noted that when the tem- 
perature is higher, potential amplitude becomes small and less efficient in trapping electrons 
because the thermal effect becomes significant. When the initial temperature is as high as 
the one corresponding to the drift velocity, the Buneman instability docs not develop. Thus, 
the attained temperature by the instability remains of the same order. The conclusion that 
surfing acceleration does not work is also valid. 

Our simulations do not treat three-dimensional effects. As long as the magnetic field 
is weak, the y- and z-dircctions are not so much discriminated, and we expect that our 
conclusions remain safe. When the magnetic field is extremely large, the situation may 
change, but at the same time the Alfven Mach number becomes small and the shock will be 
weaker. 

Our results have several implications for the electron heating and acceleration in SNR 
shocks. The attained electron temperature in simulations is comparable to the observed 
temperature of the thermal X-ray components from young SNRs. As electrons pass the front 
to downstream, they suffer further adiabatic heating, which may explain the small difference 
between the attained and observed temperatures, in addition to the Coulomb heating. This 
in turn may suggest that the degree of further heating in the upstream such as that due to 
the excitation and decay of ion acoustic modes is not significant. This point deserves further 
study. As for the acceleration of nonthermal electrons, the question remains open. Since the 
attained temperature is on the order of kilo-electron volts, it is still 4 orders of magnitude 
short of the energy necessary to enable injection into the diffusive shock acceleration. 



4. Summciry 

We performed electrostatic 2D PIC simulations to investigate electron heating and ac- 
celeration by coUisionless shocks with a high Mach number. We consider only the foot region 
by properly modeling the effects of refiected protons. We have shown that excited electro- 
static waves are oblique to the shock normal, and the electrostatic potential loses coherence 
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in the direction perpendicular to the shock normal and the magnetic field. As a result, 
surfing acceleration seen in the previous ID simulations does not occur in 2D simulations or 
in real cases. 

We are grateful to T. Tsuribe, T. Kato and Y. Pujita for discussions and suggestions 
especially for providing useful guidance and information in doing numerical simulations. 
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Fig. 1. — Time development of the spatially averaged energy density of the electric field. 
Solid, dashed, and dotted curves represent the ID simulation and the x- and y-components 
in the 2D case, respectively. 
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Fig. 4. — Energy spectrum of electrons at the final state of our simulation t=1000 u^^. 
Solid, dashed and dotted histograms represent the ID result, the 2D one with magnetic field 
and the 2D one without magnetic field, respectively. 



